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Abstract 

\& . In this work, we extend the analytic treatment of Bessel functions of large order and/or 

argument. We examine uniform asymptotic Bessel function expansions and show their accuracy 
' and range of validity. Such situations arise in a variety of applications, in particular the Fourier 

transform of the gravitational wave signal from a pulsar. The uniform expansion we consider 
O ■ here is found to be valid in the entire range of the argument. 
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1 Introduction 

The detection of gravitational waves (GW) from astrophysical sources is one of the most important 
problems in experimental gravitation today. Large laser interferometric gravitational wave detectors 
like LIGO, VIRGO, LISA, TAMA 300, GEO 600 and AIGO may in the near future open a new 
window for the study of a great variety of nonlinear curvature phenomena. 

In recent works PQI2lEj we have analyzed the Fourier transform (FT) of the Doppler shifted GW 
signal from a pulsar with the use of the plane wave expansion in spherical harmonics (PWESH), 
which has a variety of applications in many areas [UGHE]. The consequent analysis of the Fourier 
transform of the GW signal from a pulsar has a very interesting and convenient development in terms 
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of the resulting spherical Bessel, generalized hypergeometric, gamma and Legendre functions. Both 
rotational and orbital motions of the Earth and spindown of the pulsar can be considered in this 
analysis which happens to have a nice analytic representation for the GW signal in terms of the above 
special functions. The signal can then be studied as a function of a variety of different parameters 
associated with both the GW pulsar signal as well as the orbital and rotational parameters. The 
numerical analysis of this analytical expression for the signal offers a challenge for fast and high 
performance parallel computation. Recent studies of the Cosmic Microwave Background Explorer 
have raised the interesting question of the study of very large multipole moments with angular 
momentum / and its projection m going up to very large values of I ~ 1000. Such problems warrant 
an intensive analytic study supplemented by numerical and parallel computation. 

Since our FT depends on the Bessel function, a computational issue arises due to large values 
of the index or order n of the function. In the GW form of the pulsar, the Doppler shifted orbiting 
motion gives rise to Bessel functions J n ( 2 ' K f° A sin 6 ) ; where 2n f° Asmd j s large for non- negligible angle 
9. Even for sin# ~ j^ni the argument is large necessitating the consideration of large values of n. 
The motivation of this work is to extend the analysis in Watson j7] for large index, argument and 
overlapping situations. Meissel [8] has made derivations for large order Bessel functions both when 
the argument is smaller than the order and vice versa. These were studied in the GW context in an 
earlier work by us jU] . We extend this study by considering another type of expansion, also referred 
to as the "uniform" Bessel expansion, which, being an expansion in the Airy functions, allows for 
single series form, hence the name ^Hj. As an application, we will address the phenomenological 
situation of GW signal analysis of large order n (which does arise with combinations of I and m). 
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The expansion in terms of the Bessel functions has recently found a valuable application JT] 
in the identification of global parameter space correlations of a coherent matched filtering search 
for continuous GW from isolated neutron stars. The authors have done an interesting pioneering 
analysis of the global correlations and indicate that the locations of local maxima of the detection 
statistics involve Bessel functions and are dominated by the Doppler shift due to the Earth's orbital 
motion. We have made a preliminary extension of their analysis that makes use of the summation of 
infinite series involving Bessel functions and other relevant functions of the parameters and physical 
variables. Their various levels of approximation can be further improved for a comparison with exact 
numerical results. 

Captures of stellar-mass compact objects (CO) by massive black holes are important capture 
sources for the Laser Interferometer Space Antenna (LISA), the space based GW detector due to 
be launched in about a decade ^2]- Higher Harmonics of the orbital frequency of the COs arise 
in the post Newtonian (PN) capture GW model forms and contribute considerably to the total 
signal to noise (S/N) ratio of the waveform. The GW form can be decomposed into gravitational 
multipole moments which are treated in the Fourier analysis of Keplerian eccentric orbits. The 
radiation depends strongly on the orbital eccentricity e, and Bessel functions J n (ne) are a natural 
consequence of the analysis. Such calculations involving gravitational wave forms necessitate a 
very accurate evaluation of Bessel functions of large order and argument due to the importance of 
orbital eccentricity. Accurate and efficient evaluations of Bessel functions for constructing reliable 
gravitational wave templates from binary stars with arbitrary orbital eccentricity have been done 
by Pierro et al. ^3] using generalized Carlini-Meissel expansions. The methods investigated in the 
current work can also be applied to such calculations relating to LISA capture sources. 

The calculation of partial derivatives of the potential scattering phase shifts which often con- 
tain Bessel and Legendre functions of large order angular momentum /, with respect to angular 
momentum arise in a variety of scattering problems in atomic, molecular and nuclear physics. In 
particular, large values of I can arise in rainbow, glory and orbit scattering JUEI]j m black hole 
gravitational fields [121 , and in catastrophe theory 17. The analysis in our paper should help 
provide suitable approximations for large order and/or argument for the Bessel functions that arise 
in such problems. 

In section 2 we give the relevant expressions from our work on GW signal analysis to demonstrate 
the occurrence of large order Bessel functions. Section 3 presents the uniform expansion and the 
expressions associated with it. We discuss the results in Section 4 and compare them with other 
expansions and the exact values of the Bessel functions. 

2 Fourier Transform of the GW signal 

The FT for the GW Doppler shifted pulsar signal [TJ is given as follows: 
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Here a is the co-latitude detector angle and angles 8 and </> specify the direction of the pulsar 



source. Also u; = 27r/ , o; orb = (T orb = 365 days), u; rof = ^- (T rE = 1 day), B, 
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A = 1.5 x 10 11 meters is the Sun-Earth distance. 
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(Re is the radius of Earth, c is the velocity of light) and 



3 Uniform expansion of the Bessel function 

An expansion that approximates the Bessel function J v {yz) for large v, which in principle is valid 
for all real values of z, is given in JU] as 
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Ai(x) represents the Airy function. Equations IT21 and 1131 are used to define C> the choice de- 
pending on the region of interest. Equation ED is used when the argument of the Bessel function 
is less than the order, while equation is used when the argument is larger than the order. Also, 
the Ui are defined recursively by 

u (t) = 1 

u k+ i(t) = \t 2 (1 - t 2 ) u' k {t) + | / * (1 - 5t 2 ) u k {t)dt (k = 0, 1, 2, ...) 

Using the symbolic package MAPLE, we have computed these functions up to Uu, and the results 
are given below. They may, however, be calculated in principle up to any order, and MAPLE can 
perform this task fairly rapidly. 
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__ 85085 y9 _j_ 17017j.7 _ 4563 y.5 _i_ 75 f 3 
3 82944 9216 5120 1024 



_ 37182145 f 12 _ 7436429 ylO , 144001 f 8 _ 96833 ^6 _i_ 3675 f A 
4 7962624 663552 ^ 16384 40960 32768 



5391411025 ^15 i 5391411025 y!3 _ 108313205 y.11 , 250881631 y9 _ 67608983 f , 59535 f 5 
191102976 63700992 1 1179648 5898240 9175040 ~ t ~ 262144'' 



5849680962125 y.18 _ 1169936192425 y.16 i 4445922195 y!4 _ 33010308331 y!2 , 1441372804469 f 10 

6 27518828544 1528823808 4194304 47185920 6606028800 

388895895 f S i 2401245 y.6 
14680064 1 4194304 

1267709431363375 f 21 , 1774793203908725 y.19 _ 36927006432745 y.17 , 10559432785187 y!5 

7 660451885056 ' 220150628352 2717908992 905969664 1 

1602251736839 y!3 , 1007390378503 yll _ 25388505925 y9 , 57972915 ^7 
301989888 ^ 838860800 234881024 1 33554432 



_ 2562040760785380875 y24 _ 512408152157076175 y22 i 75358832548684685 y20 _ 39803268297948155 y!8 i 
8 126806761930752 5283615080448 391378894848 195689447424 ~ l ~ 

3542717254441859 y.16 _ 276439228010667 y!4 , 667955999804539 f \2 _ 928090660435 ylO , 13043905875 f 8 
28991029248 6710886400 93952409600 1879048192 2147483648 



6653619855759634132375 y27 i 1330723971151926826475 y25 _ 3128960418491082175 y23 
27390260577042432 1014454095446016 1043677052928 



35348759075759093965 y.21 _ 17618708259302571707 y.19 , 817138105244771959 y!7 _ 3739063570455884033 y.15 i 

9393093476352 6262062317568 ' 644245094400 11274289152000 ~ t ~ 

1359491937582325 y.13 _ 472414367256615 f ll , 418854310875 f 9 
30064771072 188978561024 6 17179869184 1 
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_ 4318199286388002551911375 +30 _ 4318199286388002551911375 f 2S _i_ 91891691150784941691275 y.26 
10 1314732507698036736 219122084616339456 ^ 1803473947459584 

16705838021516291703055 +24 , 19941766574064397067317 +22 _ 5227733363217471800551 f 20 , 
225434243432448 ' 300578991243264 L 139156940390400 ' 

1232816120293110459821 ^18 _ 35892416277828185849 y-16 i 45660648644355162105 +14 _ 20993386079260455 +12 
92771293593600 12884901888000 ' 148159191842816 1511828488192 

I 30241281245175 +10 
' 274877906944 1 

_ 1556514560957130010757145625 ^33 i 3424332034105686023665720375 ^31 _ 365967912305800454531456575 y.29 
11 31553580184752881664 10517860061584293888 389550372651270144 

i 201734750392525792544487385 +21 11694306169843138084657687 ^.25 i 164293183874328160710877 ^23 



+zo I iD4z«jig^a^jzgiDuriua^ +z 
129850124217090048 " 7213895789838336 ' 148434069749760 

23186185730591085896097833 y.21 , 527174389121818780771231 +.19 _ 329641577686894230674187 ^-17 i 



46756731971174400 1 3710851743744000 13469017440256000 

241770821762631191867 f 15 _ 26416375998266454375 f 13 1 1212400457192925 y.11 
107752139522048 314460325543936 2199023255552 1 

In Table [0 we compare the values generated by this expansion with the exact Bessel function 
values, as well as the other expansions (Watson's expansion, epsilon expansion and Meissel's third 
expansion) from our previous work [5] for order v = 300. (Refer to the discussion section for the 
formulae for these expansions.) We also give a comparison for a larger order, v = 50000, in Table 
121 We were not able to give exact values of the Bessel function for v = 50000 due to the immense 
computational power required, but we give values for the other expansions for comparison. 
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Table 1: Comparison of expansions for J u (x) with v = 300. 



X 


290 


295 


298 


302 


305 


310 


Watson 


0.007669 


0.027187 


0.048992 


0.084432 


0.100255 


0.057368 


Epsilon 


0.0076775 


0.027215 


0.049028 


0.084335 


0.10014 


0.057419 


Uniform 


0.007677 


0.027213 


0.27573 


-0.083138 


0.10014 


0.057419 


Exact 


0.007677 


0.027215 


0.049028 


0.084335 


0.10014 


0.057419 



Table 2: Comparison of expansions for J u (x) with v = 50000. 



X 


49800 


49900 


49999 


50001 


50100 


50200 


Watson 


3.84377 x 10~ 8 


1.02742 x 10" 4 


0.011802 


0.012444 


-0.013648 


0.001398 


Epsilon 


80.294 


0.0016699 


0.011839 


0.012444 


-0.012451 


-6.3189 


Uniform 


3.8445 x 10" 8 


0.00010275 


8.7321 x 10 i3 


3.8936 x 10 11 


-0.013648 


0.001398 


Exact 


N/A 


N/A 


N/A 


N/A 


N/A 


N/A 



The uniform expansion values were calculated using four terms in the series in equation [7| We 
note that there is a significant discrepancy between the uniform expansion values and the exact 
values close to the transition region. Temme ^HJ has noted that this is due to numerical singularities 
rather than to analytic singularities, and he has given expansions of the functions dk(C) and 
appearing in equation Using these improved formulae, we obtain the following uniform expansion 
values. 



X 


290 


295 


298 


302 


305 


310 


v = 300 


0.00767703 


0.027215 


0.049028 


0.084335 


0.100143 


0.057419 


X 


49800 


49900 


49999 


50001 


50100 


50200 


v = 50000 


3.844421 x 10~ 8 


0.00010275 


0.011839 


0.012444 


-0.013648 


0.0013977 



A similar expansion for the first derivative of the Bessel function is also given in ^U|. This 
expression is given below. 
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v k (t) = u k (t) +t(t 2 - 1) [tu'^t) + |u fc _i(t)] 

(A; =1,2, ...) 

Using MAPLE, we have calculated these functions up to V\i, although generating higher order 
expressions is trivial, as is the case with Ui. We give these results below. 
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775249765 y!4 _i_ 7176153985 y!2 _ 75861726551 yip , 
4194304 9437184 314572800 ~ t ~ 
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905969664 1 



838860800 17616076800 1 234881024 " 33554432" 
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126806761930752 587068342272 1174136684544 195689447424 

!57093219318079 y!6 i 296916207863309 y!4 _ 29041565208893 y!2 , 146540630595 ylO _ 14783093325 yg 
9663676416 ' 6710886400 3758096384 1 ' 268435456 2147483648 

6904699850316601458125 y27 _ 461679745093525633675 y25 i 29412227933816172445 y23 
' 27390260577042432 338151365148672 ' 9393093476352 

37073088786771732695 y21 , 6190356955971173843 y!9 _ 519996976064854883 y.17 _i_ 3996930023590772587 y!5 
9393093476352 2087354105856 1 386547056640 1 11274289152000 L 

1468251292588911 y!3 , 517406211757245 yll _ 468131288625 y9 
30064771072 188978561024 17179869184 

_ 4464578923214714502823625 y30 i 1491741571661309972478475 y28 _ 286485860646564818213975 y26 i 
1314732507698036736 73040694872113152 5410421842378752 ~ l ~ 

17416724745836133903185 y24 _ 2318810066751674077595 y22 i 16487466760916641832507 y20 
225434243432448 33397665693696 417470821171200 

45614196450845087013377 y.18 , 12736018679229356269 y.16 _ 49042918914307396335 y.14 i 22818897912239625 y!2 
3246995275776000 ^ 4294967296000 148159191842816 ~ l ~ 1511828488192 

_ 33424574007825 ylO 



33424574007825 f 
274877906944 

v n = 



274877906944 

1604407316678887857241980875 y33 _ 392956135061308232223935125 y31 i 1136426675054854043018733575 y29 
31553580184752881664 1168651117953810432 1 *~ 1168651117953810432 

209347382482809784715977475 y27 , 4057208263006803008962871 y25 _ 7721779642093423553411219 y23 , 
129850124217090048 * 2404631929946112 1 6679533138739200 ~ l ~ 

24317219180863821793468459 y21 _ 185223434015774166216919 y!9 _i_ 109880525895631410224729 .17 
46756731971174400 1236950581248000 ~ t ~ 4233119766937600 

258444671539364377513 y!5 _i_ 28529686078127770725 y!3 _ 1327867167401775 yll 
107752139522048 314460325543936 2199023255552 1 

In TableElbelow, we compare the values of J' v (x) obtained from the expansion with the exact values, 
for v = 300. A similar table is not given for v = 50000 since we are not able to provide exact values 
for comparison, due to the large value of the index. 
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Table 3: Comparison of uniform expansion of J' u (x) with exact values for is = 300 



X 


290 


295 


298 


302 


305 


310 


Uniform 


0.0021909 


0.0059599 


-0.010098 


-0.019458 


0.0014413 


-0.018938 


Exact 


0.0021909 


0.0059597 


0.0084469 


0.008023 


0.0014411 


-0.018938 



4 Discussion 



We find that the uniform expansions are accurate to a higher domain of validity than the expansions 
considered in our earlier work [Oj, thereby increasing their utility especially in the problematic z ~ 1 
or "transition" region. We compare the uniform expansion with the exact Bessel function values 
in the transition region in Figures El and H3 at the end of this section, which are plotted for order 
is = 300. It is seen that the uniform expansion agrees with the exact Bessel function curve for the 
entire range of the argument. 

Previously, we presented two methods that are geared to work specifically in the transition 
region 0. First, we present the results by Watson [7j. For the case of the argument being less than 
the order, he obtained, through the use of contour integration, 



T . , tanh a 
J v [ysech\a,)) = j=- exp 



7T 



V3 



( tanh a H — tanh 3 a — a ) 

v 3 ; 



Ki 

3 



-is tanh 3 a 
3 



+ 361V 1 exp[^(tanhc>! — a)] 



(17) 



where \9i\ < 1. Similarly, for the case when the argument is greater than the order, he derived the 
following: 



(18) 



J u {is sec (3) = | tan /3 cos [is (tan/3 — | tan 3 (3 — (3)] (j-i + J± 

+3-5 tan p sin [u (tan (3- § tan 3 (3-0)] (j_i-Ji) +2A6 2 v~ l 



where \9^\ < 1 and the argument for the Bessel functions J.i is |z/tan 3 /3. The great advantage 

3 ^ 

of these formulae is that they have error bounds given. We have found that \9\\ and \9 2 \ do not 
generally have constant values for different values of the argument, but they do nevertheless stay 
below 1. We have plotted the values of 9\ and # 2 such that Watson's expansion matches the exact 
value, for a range of arguments with v = 300 (Figure Q). From these graphs, it is clear that \9i\ 
stays below ~ 0.0052 and \9 2 \ stays below ~ 0.0016. 

On the other hand, Debye ^H] introduced what we will term as the "e expansion" . The idea is 
motivated by introducing a small parameter e, such that is = x(l — e), where v denotes the order 
and x is the argument of the Bessel function. This expansion is given below: 



J v (x) ~ — B m (ex) sin 



m=0 



(m 



lk 



a x ym+l)/3 



(19) 



We extended this analysis of Debye by 5 orders to the terms B m (ex),m = 0,1, 2,. .15 0. It 
should be remarked that the single value of the Bessel function when the argument equals the order 
(x = is) was also considered in this earlier study which extended the previous work of 20J by two 
orders. We observed that inclusion of the higher order terms leads to highly accurate values, and 
so the singular case, for the uniform expansions (when x = is) is well under control. 
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Figure 1: Plots of calculated values of (a) Q\ and (b) 62 vs. the argument of the Bessel function, 
for order v = 300. 

Now, to illustrate the applicability of these three methods to the transition region (z ~ 1), we 
present Figures 0] and which are plotted for the problematic regions relevant for GW phenomenol- 
ogy (when the order is as large as v = 1000000). Watson's expansion and the e expansion both 
show remarkable ability in capturing the functions in the transition region. 

We were not able to make exact comparison for such a large order, obviously due to massive 
computer times required. However, in Figures 0] and 03 we observe strong evidence that the proposed 
uniform Bessel asymptotic expansions are appropriate for GW signal analysis. Here, we note the 
transition region starting at values of the argument at 1,000,000 and going up to 1,000,200. In 
this region, the e expansion, Watson's formula and the uniform expansion closely coincide with one 
another. As we found earlier [H| the e expansion breaks down earlier as the argument leaves the 
transition region; however, all three methods coincide in a certain region indicating that we have 
consistent methods that work for values relevant to GW analysis. The uniform Bessel expansion 
is fairly easy to implement computationally and indicates good stability for rather large values of 
the argument. To illustrate the type of values a GW pulsar FT would require, we present Figure 
El Here, we plot the uniform Bessel expansion for values ranging from 1,000,200 to 32,500,000, 
which are relevant for GW phenomenology This appears as a black band and is a continuous 
function which indicates oscillations tightly bunched together, and happens to be identical to a 
graph obtained earlier by using Meissel's second expansion jH] . It is noteworthy that the method is 
stable and shows consistent behaviour over an extreme range of values for the argument. 

The zeroes (roots) of the Bessel function will help to identify the sky locations in terms of 9 
where the signal strength is zero. Zeroes of Bessel functions arise in a variety of diffraction problems. 
Derivatives of Bessel functions are also important in a multitude of applications and are evaluated 
using the uniform expansions. 

It is worth mentioning here a procedure developed by Baruth j^] which can evaluate Bessel 
functions for large index and argument, for example in the range of 5 million. The procedure per- 
forms evaluations in less than 10 ms, and the author believes that it can be improved even further. 
As an example, he has calculated J1000000 (999995) = 4.267834146855037 x 10~ 3 , which is close to 
our value of 4.267834146855055 x 10~ 3 , an agreement of 13 decimal places. 

Number-theoretic study of the denominators of Ui 
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We note here an interesting pattern regarding the functions Ui used in the uniform expansion. 
Using MAPLE we undertook a study of the denominators occurring in the Ui mainly because their 
prime number factorizations appear to have some surprising properties. 

Let di denote the common denominator of the terms of Ui for i = 0, 1, . . . and consider the prime 
number factorizations of the first eleven d^s : 

i di 



1 

1 2 3 3 

2 2 7 3 2 

3 2 10 3 4 5 

4 2 15 3 5 5 

5 2 18 3 6 5 7 

6 2 22 3 s 5 2 7 

7 2 25 3 9 5 2 7 

8 2 31 3 10 5 2 7 

9 2 34 3 13 5 3 711 

10 2 38 3 14 5 3 7 2 11 

Powers of 2 

Let cij denote the exponent of 2 in the prime number factorizations of di for i — 0, 1, Then the 

following formula for seems to be valid 

at = log 2 (21) 

where gi is Gould's sequence, see |22], defined as the highest power of 2 dividing the central binomial 
coefficient ( 2l ). All terms of Gould's sequence gi are powers of two and therefore a« will be an 
integer. This remarkable identification was done via N. J. A. Sloane's On-Line Encyclopedia of 
Integer Sequences, see |2Sj. 

The table below lists the first eleven values given by formula (|21|) which agree with the values 
given in the previous table (J2UJ) . 



i 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


9i 


2° 


2 1 


2 1 


2 2 


2 1 


2 2 


2 2 


2 3 


2 1 


2 2 


2 2 


di 





3 


7 


10 


15 


18 


22 


25 


31 


34 


38 



Powers of 3 

Let bi denote the exponent of 3 in the prime number factorizations of di for i — 0, 1, . . .. Then the 
following formula for bi seems to be valid 

bi = i + exponent of the highest power of 3 dividing i\ (22) 

The table below lists the first eleven values given by formula (|22j) which agree with the values 
given in (J2Dj). 



% 





1 


2 


3 


4 


5 


6 7 


8 


9 


10 


exponent 











1 


1 


1 


2 2 


2 


4 


4 


bi 





1 


2 


4 


5 


6 


8 9 


10 


13 


14 



10 



(20) 



These number theoretic properties could potentially be useful toward finding an analytic solu- 
tion to the integro-differential equation describing the Bessel function for the uniform asymptotic 
expansions. 
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Figure 2: Uniform Bessel expansion and actual 
Bessel function graphed for argument x and 
order v = 300 near the transition region. Solid 
line indicates actual Bessel function values and 
circles indicate values given by the expansion. 




Figure 3: Uniform Bessel expansion and actual 
Bessel function graphed for argument x and 
order v = 300 near the transition region. Solid 
line indicates actual Bessel function values and 
circles indicate values given by the expansion. 



9.99759+05 9.998e+05 9.9985e+05 9.999e+05 9.9995e+05 



Figure 4: Comparison of uniform expansion, e 
expansion and Watson's formulae for argument 
x < 1, 000, 000 and order v = 1, 000, 000. Solid 
line indicates uniform Bessel expansion values, 
diamonds represent e expansion and circles in- 
dicate values given by Watson's formula. 




Figure 5: Comparison of uniform expansion, e 
expansion and Watson's formulae for argument 
x > 1, 000, 000 and order v = 1, 000, 000. Solid 
line indicates uniform Bessel expansion values, 
diamonds represent e expansion and circles in- 
dicate values given by Watson's formula. 
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Figure 6: Plot of uniform Bessel expansion for argument x ranging from 1,000,200 to 32,500,000 for 
order 1,000,000. 



5 Conclusions 

In this work we examined, using symbolic computation, the accuracy and the range of validity of 
uniform asymptotic expansions of Bessel functions relevant for GW signal from pulsars and found 
them to be valid over the entire domain of the argument, including the problematic "transition 
region" . We compared these expansions with the fractional transition Bessel and the Debye expan- 
sions in such regions and found all three expansions to be applicable in the transition region. It was 
shown that there is a region where all three expansions agree, which is indicative of the accuracy 
of the uniform expansion in this domain. 

The functions used by Bessel in 1824 in connection with planetary motion have found innu- 
merable applications as earlier studies by Bernoulli, Lagrange, Carlini, Laplace, Poisson, Euler and 
many others [2j have shown. Further problems will also bring the ubiquitous Bessel functions to 
more diverse applications. Such applications may require accurate calculations of the functions for 
large order and argument, and asymptotic analysis will also be relevant for other special functions. 

6 Acknowledgements 

We are deeply grateful to SHARCNET (Shared Hierarchical Academic Research Computing Net- 
work) and NSERC for grant support. We are indebted to Dr. Nico Temme (CWI, Amsterdam) 
for suggesting the importance of this study and for his invaluable advice on the evaluation of the 
uniform expansion, and also to Dr. Dan Baruth for bringing to our attention his work on evaluation 
of Bessel functions of large order and argument. 




13 



References 

[1] K. Jotania, S.R. Valluri and S.V. Dhurandhar, Astron.Astrophys. 306:317,1996. 

[2] S.R. Valluri, F.A. Chishtie, et al. Class.Quant.Grav.l9:1327-1334,2002, Erratum-ibid.l9:4227- 
4228,2002. 

[3] Valluri S. R., Chishtie F. A. and Vajda A. Classical and Quantum Gravity, 23, 3323-3332 
(2006). 

[4] MacPhie, R.H. and Wu, K.L., IEEE Trans, on Antennas and Propagation, vol.51, No. 10, 
pp.2801-2805, 2003. 

[5] Thorne, K.S., "Multipole expansions of gravitational radiation", Reviews of Modern Physics, 
vol.52, pp. 299-339, 1980. 

[6] Allen, B. and Ottewill, A.C. Detection of Anisotropies in the Gravitational- Wave Stochastic 
Background. Phys. Rev. D 56:54563 

[7] Watson, G.N., A treatise on the theory of Bessel Functions, Cambridge University Press, 1958; 
Proc. Camb. Phil. Soc. vol. XIX, 96, 1918. 

[8] Meissel, D.F.E,"Neue Entwickelungen ueber die Bessel'schen Functionen" ylstr. Nach. vol. 
CXXIX col. 281-284, 1892; Astr. Nach. vol. CXXX , "Weitere Entwickelungen ueber die 
Bessel'schen Functionen", col. 363-368, 1892; "Einige Entwickelungen die Bessel'schen I- 
Functionen betreffend" vol. CXXVII, col. 359-362, 1891; "Beitrag zur theorie der allgemeinen 
Bessel'schen Functionen" vol. CXXVIII (1891), col. 145-154. 

[9] Chishtie, F.A., Valluri S.R., Rao K.M., Sikorski, D.J., and Williams, T., IEEE Computer 
Society Press, Proc. 19th Annual Symposium, High Performance Computing Systems & Ap- 
plications. (HPCS 2005), Guelph, Ont., Canada. 

[10] Abramowitz, M. and Stegun, LA., Handbook of Mathematical Functions, National Bureau of 
Standards, 1964. 

[11] Prix, R. and Itoh, Y. "Global parameter-space correlations of coherent searches for continuous 
gravitational waves". Classical Quantum Gravity 22:S1003-S1012 (2005). 

[12] Peters P. C. and Mathews J., "Gravitational Radiation from Point Masses in a Keplerian 
Orbit", Phys. Rev. vol. 131, pp. 435-439, 1963; Barack, L. and Cutler, C, "LISA Capture 
Sources: Approximate Waveforms, Signal-to-Noise Ratios, and Parameter Estimation Accu- 
racy" Phys. Rev. D vol. 69, 082005, 2004. 

[13] Pierro, V., Pinto, I.M., Spallicci, A.D., Laserra, E. and Recano, F., Mon. Not. R. Astron. Soc. 
325:358-372 (2001). 

[14] Berry, M.V. and Mount, K.E., Rep. Progr. Phys. 35, 315 (1972). 

[15] Ford, K.W. and Wheeler, J.A., Ann. Phys. 7:259, 287 (1959); Romo, W.J. and Valluri, S.R., 
Physica Scripta 71:1-10 (2005). 

[16] Regge, T. and Wheeler, J.A., Phys. Rev. 108, 1063 (1957). 



14 



[17] Poston, T. and Stewart, I. Catastrophe theory and its applications, Pitman Publishing, London, 
1978. 

[18] Temme, N.M., Numerical algorithms for uniform Airy-type asymptotic expansions, Numerical 
Algorithms 15, 207 - 225 (1997). 

[19] Debye, P., "Naeherungsformeln fuer die Zylinderfunktionen fuer grosse Werte des Arguments 
und unbeschraenkt veraenderliche Werte des Index", Math. Ann., vol. LXVII pp. 535-558, 1909. 

[20] Airy, J.R. "Bessel and Neumann Functions of Equal Order and Argument" Phil. Mag. (6) vol. 
XXXI, 520-528, 1916. 

[21] Baruth, D., Bessel Functions of Large Orders and Arguments. World Wide Web URL 
www. iging.com/transcendental/ J_ 1000000 .htm 

[22] J. -P. Allouche and J. Shallit, The ring of k-regular sequences, Theoretical Computer Sci., 98 
(1992), 163-197. 

[23] N. J. A. Sloane, The On-Line Encyclopedia of Integer Sequences. World Wide Web URL 
www . research . att . com/ nj as/sequences/ 



15 



